Collective structural evolution in neutron-rich Yb, Hf, W, Os and Pt isotopes 
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An interacting boson model Hamiltonian determined from Hartree-Fock-Bogoliubov calculations 
with the new microscopic Gogny energy density functional DIM, is applied to the spectroscopic 
analysis of neutron-rich Yb, Hf, W, Os and Pt isotopes with mass A ~ 180 — 200. Excitation 
energies and transition rates for the relevant low-lying quadrupole collective states are calculated 
by this method. Transitions from prolate to oblate ground-state shapes are analyzed as a function 
of neutron number N in a given isotopic chain by calculating excitation energies, S(E2) ratios, 
and correlation energies in the ground state. It is shown that such transitions tend to occur more 
rapidly for the isotopes with lower proton number Z, when departing from the proton shell closure 
Z = 82. The triaxial degrees of freedom turn out to play an important role in describing the 
considered mass region. Predicted low-lying spectra for the neutron-rich exotic Hf and Yb isotopes 
are presented. The approximations used in the model and the possibilities to refine its predictive 
power are addressed. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60.Jz 
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I. INTRODUCTION 

The study of the origin of nuclear deformation and 
its evolution as a function of proton and neutron num- 
bers has attracted considerable theoretical interest from 
a large variety of viewpoints [H-[T3j . Experimentally, low- 
lying spectroscopy provides a very powerful source of in- 
formation that allows one to establish signatures correlat- 
the nuclear shape evolution with the energy spectra 



Among many other nuclear structure models, self- 
consistent mean-field methods, based on microscopic en- 
ergy density functionals (EDFs), have provided both 
accurate and universal descriptions of different nuclear 
intrinsic properties including binding energies, ground- 
state deformations, density distributions, low-lying one- 
quasiparticle configurations, as well as the way nuclear 
shapes evolve with the number of nucleons @, H, @- 
Inl 12414 2 3 . Popular EDFs are the non-relativistic Skyrme 
[6[|29|,|3(| and Gogny 3ll32| ones, as well as relativistic 
mean- field Lagrangians 33]. To describe nuclear spec- 
troscopy one should go beyond the mean-field approxi- 
mation to take into account the restorations of broken 
symmetries and/or the configuration mixing of intrinsic 
states in the spirit of the generator coordinate method 

(GCM) [HUMS- In thiS kind ° f studies calculations 
may become computationally much more demanding and 

time consuming than the underlying mean field, particu- 
larly when triaxial degrees of freedom are included in the 
analysis. 

A sound approximation to the full GCM configura- 
tion mixing and/or the symmetry restoration is the five- 
dimensional collective Hamiltonian with quadrupole de- 
grees of freedom where both rotational and vibrational 



mass parameters are determined from the constrained, 
self-consistent mean-field calculations with a given EDF 
and the collective potential is derived by the zero-point 
energy correction to the total mean- field energy (e.g., 

Alternatively, nuclear dynamics and spectroscopic 
quantities can be approximated by introducing appro- 
priate bosonic degrees of freedom. The interacting bo- 
son model (IBM) [42j can be regarded as a nice exam- 
ple for this, and has been exploited in a large number 
of phenomenological studies focusing on the low-lying 
spectrum of medium- heavy and heavy nuclei 42]. The 
simplest version of the IBM is built on monopole s and 
quadrupole d bosons, which reflect the collective J 71 ' = + 
and 2 + pairs of valence-shell configurations, respectively 
43]. Nevertheless, since the IBM itself should have a cer- 
tain microscopic foundation, a Hamiltonian of IBM has 
been derived conventionally from the shell-model config- 
uration |43| , and more recently from EDF-based calcula- 
tions [44| • These mapping methods have been applied to 
realistic cases involving a variety of situations covering 
from nuclei with modest quadrupole deformation includ- 
ing 7 -unstable ones [43-47], to strongly deformed rota- 
tional nuclei (48l - l53j . Also quantum-mechanical correla- 
tion effects in the ground state have been considered (4(| . 
Starting from the constrained Hartree-Fock-Bogoliubov 
(HFB) theory with the D1S 54j parametrization of the 
Gogny functional, the method of 44J was used for the 
spectroscopic analysis of Pt isotopes ^7[, and some of 
Os and W isotopes [53| . 

In this paper we apply the mapping procedure of [44| 
to the mass regi on A ~ 180 — 200, extending the analysis 
made in [U [53| to the neighboring exotic Hf and Yb nu- 
clei. An additional motivation is to explore some possibil- 
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ities to refine the predictive power of the method for the 
considered mass region. Although the D1S parametriza- 
tion of the Gogny force is considered as a global EDF able 
to describe many low-energy nuclear data with reason- 
able predictive power (see, for example, Refs. (4ll.l55l.l56j 
and references therein), we have preferred to use in this 
paper the Gogny-DIM functional [5^ |. Systematic ex- 
plorations of different nuclear phenomena [55l . in- 
cluding properties of odd nuclei computed within the 
equal- filling approximation [24Tj26| , suggest that the new 
parametrization of the Gogny-EDF is as good as the stan- 
dard D1S, a fact that we intent to confirm in this paper. 

It is well known that the nuclei in the mass region 
A ~ 180 — 200 exhibit a transition between prolate and 
oblate equilibrium shapes as a function of the nucleon 
number, with the critical point around N w 116 having 
a pronounced 7 softness l5a - l61j |. These facts make the 
region a potential testing ground to understand the de- 
formation properties of atomic nuclei. The evolution of 
the nuclear ground states in this mass region has been 
investigated recently with the constrained self-consistent 
mean- field method with microscopic EDFs [9|, |56|, l62j . 
Both the (constrained) Hartree-Fock+BCS (HF+BCS) 
and the HFB approximations have been used to com- 
pute energy surfaces with quadrupole degrees of freedom 
in order to give a microscopic insight into shape transi- 
tions [j| [55l |56| . It was shown in these studies that the 
triaxiality is an important ingredient to describe the evo- 
lution from prolate to oblate shapes, irrespective of the 
types of the EDFs used. 

It should be kept in mind that Pt, Hg, and Pb iso- 
topes are well known [§3[ for the spectacular coexistence 
of different low-lying configurations based on different in- 
trinsic deformations as observed in their low lying spec- 
trum. There are a number of works aimed at the un- 
derstanding of the shape coexistence phenomenon in this 
region in terms of both EDF-based microscopic calcula- 
tions [H, H3, Hl| and phenomenological models [66T-l68| . 

The paper is organized as follows: In Sec. [Til a short 
outline of the theoretical framework is given. Section Hm 
presents the energy surfaces, ground-state correlation en- 
ergies, moments of inertia for the rotational bands, low- 
lying spectra, and the B(E2) systematics for the consid- 
ered isotopes chains. Section |IV] is devoted to the con- 
cluding remarks and work perspectives. 



II. THEORETICAL PROCEDURE 

The analysis starts with a constrained HFB calculation 
using the Gogny-DIM EDF. The constraints in this case 
refer to the mass quadrupole moments which are associ- 
ated with the quadrupole deformation parameters /3 and 
7 in the geometrical model [l|. The set of constrained 
HFB calculations, for each collective coordinate (/3/y), 
provides the total HFB energy (denoted by Eufb(P : l))- 
For calculation details the reader is referred to [43, [56[ . 

In other studies solving the five-dimensional collective 



Hamiltonian [39h41] . the collective potential energy sur- 
face is obtained by subtracting the zero-point energies 
for both rotational and vibrational motions from the con- 
strained HFB energy surface. This corrected energy sur- 
face should be viewed as a collective potential energy sur- 
face. In the present work, the constrained HFB energy 
surface and the corresponding boson energy surface are 
compared, and they will be referred to simply as energy 
surface. Note that, as the total energy is considered, all 
ingredients including those relevant to kinetic terms are 
supposed to be taken into account to a good extent. 

Each point of the Gogny-HFB energy surface 
£"hfb(/3,7) is mapped onto the corresponding point on 
the bosonic energy surface, denoted by -Eibm(/?b, 7b) 
with Pb and jb being the deformation parameters for 
the boson system, in such a way that the bosonic energy 
surface fits the fermionic one [441 ]. In this paper we con- 
sider the proton-neutron interacting boson model (IBM- 
2) (43[ because it reflects better the microscopic picture 
than the original version of the IBM without distinction 
of the proton and the neutron degrees of freedom (of- 
ten called IBM-1). In what follows we denote the IBM-2 
simply as the IBM, unless otherwise specified. The IBM 
energy surface is obtained as the expectation value of 
a given boson Hamiltonian 69] in terms of the coher- 
ent state |$(/3b, 7b)). The coherent state represents the 
intrinsic wave function of the boson system, and is char- 
acterized by the deformation variables Pb and 73. In 
principle, proton and neutron bosons might have differ- 
ent values of the deformation parameters, but since pro- 
ton and neutron systems are supposed to attract each 
other strongly in medium-heavy and heavy deformed nu- 
clei, the deformations of proton and neutron systems can 
be taken the same to a good approximation. 

If the separability of the mapping along the /3 and the 
7 directions is assumed, one can consider the relation be- 
tween the IBM and the geometrical deformation variables 
|46[. It was shown [69| that, in general terms, the 
bosonic and the geometrical /3s are proportional to each 
other and that the proportionality coefficient coincides 
with the ratio of the total nucleon number to the valence 
nucleon number counted from the nearest closed shells. 
We exploit this relation and assume that /3b = C^/3, with 
Op being a numerical coefficient 44]. The typical range 
of the C/3 value turns out to be approximately 5 ~ 10, 
which is about the same order of magnitude as the actual 
ratios of the total nucleon number to the valence nucleon 
number. Regarding the triaxial parameter 7, the identi- 
fication 7b — 7 seems valid as indeed both geometrical 
and IBM 7's have the same meaning, ranging from to 
60 degrees. 

We adopt the IBM Hamiltonian of the following form: 

Hjbm = thd + kQk ■ Q v + ah ■ L, (1) 

where the first term fid = h^+ntiv with hd p — d p -d p (p — 
7r or v) is identified as the <i-boson number operator. The 
second term on the right-hand side of Eq. ([TJ stands for 
the quadrupole-quadrupole interaction between proton 
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and neutron systems, with Q p = s fl d p + dps p + Xp[d p d P Y 2 ^ 
being the quadrupole operator for proton or neutron sys- 
tems. The third term (denoted by LL term, hereafter) is 
relevant to the moment of inertia of the rotational band. 
L = L K + L v is the angular momentum operator for the 
boson system with L p = VlOldpdp]^. 

The form of the Hamiltonian Hibm in Eq. (p} is not the 
most general, but embodies all essential features of the 
low-lying quadrupole collective states. A more general 
IBM Hamiltonian with up to two-body interactions con- 
tains many more terms than those considered here. How- 
ever, these additional terms are supposed to be of little 
importance, and their implementation would increase the 
number of parameters, which makes the problem quite 
complicated. 



TABLE I: The parameters for the IBM Hamiltonian Hibm 
of Eq. (TTJ, as well as the coefficient Cp, obtained from the 
mapping of HFB to IBM energy surfaces for the considered 
Yb, Hf, W, Os and Pt nuclei with N = 110 — 122. 
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The parameters contained in the first two terms of the 
Hamiltonian Hibm in Eq. (p}, e, k, x-n and Xv, as well 
as the coefficient Cp, are fixed using the fitting method 
of Ref . [46| . The LL term contributes to the energy sur- 
face in the same way as the d -boson number operator, 
but with a different coefficient, 6a. Hence, the a coeffi- 
cient cannot be fixed only by the mapping of the energy 
surface. A further step is required, in order to incorpo- 
rate specific non-zero angular frequency features of the 
rotational cranking. The a value is determined by the 
procedure of Ref. [52[ , where the cranking moment of in- 
ertia was compared between fermion and boson systems. 

We then calculate the moment of inertia for the 2\ 
excited state by the Thouless and Valatin (TV) formula 
M, 



Jtv = 3/ E 1 



(2) 



Here, E~ ( stands for the 2{ excitation energy obtained 
from the self-con sistent cra nking method with the con- 
straint (J x ) = \J L(L + 1), where J x represents the x- 
component of the (fermion) angular momentum operator 
[56| . In [111 , the Inglis-Belyaev formula [TJ [ill turned 
out to be valid for the rotational regime, but the present 
TV moment of inertia appears to be more general. 

For the boson system, we calculate the moment of in- 
ertia of the intrinsic (coherent) state, denoted by Jibu, 
using the cranking formula of Ref. 



Jibm(Pb,1b) = lim — 

ui— >0 UJ 



i w BllB )\L x \m3 B , 1B )) 

Wb>1b)Wb,1b)) 



(3) 



where w and L x stand for the cranking frequency and the 
x-component of the boson angular momentum operator, 
respectively. 

While Jibm has six parameters e, k, x-k, Xv, and a, 
all of them but a are already fixed by the energy-surface 
analysis. The a value for each nucleus is obtained so 
that the Jibm value at the equilibrium point, where the 
boson energy surface E\bm{Pb, 7b) is minimal, becomes 
identical to the J7tv value at its corresponding energy 
minimum. 

The values of all derived IBM parameters are summa- 
rized in Table HI When diagonalizing the Hamiltonian in 
Eq. ([1]), the e parameter is shifted by Ae = 6a. The e 
value listed in Table U is the one with this shift. 

The diagonalization of the IBM Hamiltonian, which is 
parametrized by the set of interaction strengths summa- 
rized in Table HI generates the energies and the wave 
functions of the excited states. Diagonalization is per- 
formed in the boson M -scheme basis, where M denotes 
the z-component of the boson angular momentum oper- 
ator. With the eigenvectors of the Hamiltonian H\bm, 
the £?(E2) value is calculated: 



1 



2L + 1 



\(L'\\f^\\L)f 



(4) 



where L and L' are the angular momenta for the initial 
and the final states, respectively. In the present work 
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the E2 operator is given as T^ E2 ' = e^Q* + e v Q v , where 
Q p coincides with the quadrupole operator in Eq. ([1}, 
and thus the same values of the Xt an d Xv parame- 
ters as those listed in Table U are used in calculating 
the B(E2) values (so-called consistent-Q formalism (cf. 
(4^|)). The boson effective charges for protons and neu- 
trons are taken the same, namely e„ = e v . 



III. RESULTS AND DISCUSSION 

A. Energy Surfaces 

Figure [1] shows the mapped IBM energy surfaces for 
Yb, Hf, W, Os and Pt isotopes with 112 < N < 120. 
Each energy surface is plotted in terms of /3(= Ps/Cp) 
and 7(= 7b) up to 2 MeV from its absolute minimum, 
since most of the quadrupole collective states are within 
this range. Note that the IBM energy surfaces for N 
1 10 and 122 are not drawn as they are similar to those 
for N — 112 and 120 nuclei, respectively. The Gogny- 
D1M energy surfaces are not shown as they do not differ 
substantially from the ones depicted in [561 ] with Gogny- 
D1S. 

For all the isotopes but the Pt ones, the energy min- 
imum shifts from the prolate (7 = 0°) to the oblate 
(7 = 60°) sides as the number of neutrons increases, pass- 
ing through the most notable 7-soft nuclei with N w 116. 
The derived x-n an d Xv values for many N — 116 isotones 
then satisfy Xn + Xv ~ 0, as summarized in Table HI This 
choice of the x parameters is at the origin of the almost 
totally flat topology of the energy surface in the IBM- 
2, as seen for example in 192 0s nucleus in Fig. [TJ The 
change in the topology of the energy surface is an evi- 
dence of prolate-to-oblate shape/phase transition, which 
becomes sharper for smaller Z. The Gogny-DlS energy 
surfaces reported in [U HH were somewhat steeper in 
both /3 and 7 directions than the present Gogny-DIM 
ones. 

A difference is apparent between the energy surfaces 
of the Pt isotopes and those of the others. For the Pt 
isotopes, the variation of the energy surface takes place 
much moderately. Such slow structural transition in Pt 
isotopes was also observed in the case of the D1S func- 
tional [U H|| . While a certain quantitative difference is 
observed between the two Gogny functional results, the 
conclusion does not change. 

It should be noted that the Gogny-HFB calculation 
suggested shallow triaxial wells for the transitional, N = 
116 Os and W nuclei [56j]. In contrast, the mapped IBM 
energy surfaces in Fig. [TJ are flat in the 7 direction, as 
the only 7-dependent term of the bosonic energy surface 
is proportional to cos 37. This is the case as long as the 
boson Hamiltonian contains up to two-body interactions. 
Only when a three-body (so-called cubic) term is consid- 
ered, a stable minimum at a 7 value different from 7 = 
and 60 degrees is obtained (7J, [75| . 



B. Correlation energies 

We next discuss a signature for a shape transition from 
a simple perspective. To do this we consider the following 
quantity that will be called correlation energy hereafter, 
which was already introduced in Ref. (46|: 

^Corr = £obm(0i ) — (-fflBM)min, (5) 

where the first term £'ibm(0^") is the eigenenergy of the 
IBM Hamiltonian, Eq. (TJJ), for the L 7 ' = + ground state, 
and the second term (-ffiBM)min denotes the minimum 
value of the IBM energy surface, that is obtained by the 
variation with respect to /3 and 7. 

In the self-consistent mean-field calculation with a 
given EDF (e.g., Ref. [Hj]), the quantum-mechanical ef- 
fect can be extracted by comparing the minimum value 
of the total energy surface of the mean field with the 
17 = + eigenenergy resulting from the restoration of 
the broken symmetries and the configuration mixing. 
For calculations of correlation energies by mapping the 
EDF theory into shell model like interactions, includ- 
ing quadrupole and pairing correlations, the reader is 
referred to |76j |. 

In the present study, all correlation effects can be in- 
cluded by the diagonalization of the boson Hamiltonian, 
and the energies and the wave functions of the states 
with good angular momentum and particle number can 
be generated. Thus, the quantity defined in Eq. (O con- 
tains correlation energies coming from symmetry restora- 
tion and configuration mixing and is similar to the equiv- 
alent quantity discussed in GCM studies. 

The behavior of Eq oxx with neutron number correlates 
well with the underlying shape transition. Figure [5] shows 
that for each considered isotopic chain the correlation 
energy is maximal in magnitude at the neutron number 
N ~ 116, which corresponds to the transition point of 
the prolate-to-oblate shape transition, and decreases as 
the neutron shell closure N = 126 is approached. This is 
consistent with the overall systematic trend of the under- 
lying energy surface in Fig. [TJ These features have been 
recognized in the GCM studies (e.g., in 38]) also. For 
the Pt isotopes, the magnitude of -Ecorr decreases with 
N, indicating that a clear transition is not expected for 
these nuclei. 

Compared with the analysis by the GCM configuration 
mixing using e.g., a Skyrme functional [34| for the same 
mass region as considered here, the magnitude of the 
present correlation energy -Ecorr is rather small, whereas 
the qualitative features mentioned above do not contra- 
dict the GCM results. 

In comparison to some rare-earth nuclei such as Nd- 
Sm-Gd isotopes, where a distinct first-order shape tran- 
sition is observed [l3| , the shape transition occurs rather 
moderately in the considered mass region. Thus, con- 
trary to -Ecorr in Fig. [21 any drastic change with nucleon 
number is not expected in some other quantities in the 
ground state, like two-nucleon separation energies. 
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FIG. 1: (Color online) The IBM energy surfaces for the considered Yb, Hf, W, Os and Pt isotopes with N = 112 — 120, obtained 
by the mapping from the Gogny-DIM energy surface, depicted within 0^/3^ 0.4 and 0° ^ 7 ^ 60° up to 2 MeV excitation 
from the minimum. Contour spacing is 100 keV. 



C. Moments of inertia 

Based on the analysis in Sec. IIII B[ we discuss to what 
extent the moment of inertia is affected by the configu- 
ration mixing due to the diagonalization of IBM Hamil- 
tonian. The effect is most nicely illustrated in the W 
isotopes, for which relatively many experimental spec- 
troscopic data are available. 

We show in Fig. [3] the moments of inertia of W iso- 
topes, calculated by the cranking formula for the coher- 
ent state Jibm in Eq. ([3]) and those taken from the 2^ 
eigenenergies of the IBM and the experimental 2^~ ex- 
citation energies [77| using the rotor formula L(L + 1). 
Note that the cranking moment of inertia of the IBM is, 
due to the correction by the LL term, set identical to the 
TV moment of inertia. Thus the TV moment of inertia 
is not depicted in Fig. [3] 

The experimental moment of inertia decreases with 
N and the slope of this decrease appears to change at 
N = 116. This change suggests a gradual shape transi- 
tion. The moment of inertia of the IBM intrinsic state, 
in contrast, decreases smoothly with the exception of the 



kink at N = 114. Perhaps such a kink reflects a detailed 
shell structure irrelevant to the present work. However, 
the kink is eliminated in the moment of inertia after di- 
agonalization, which falls on the same systematics as the 
experimental data. 

It appears that, from Fig. [3J the cranking moment of 
inertia still works for the nuclei N =110 and 112, for 
which one cannot see any difference from the moment of 
inertia taken from the IBM eigenenergies. In the tran- 
sitional region of 114 ^ TV ^ 118, where according to 
Fig. [2] a large amount of correlation energy should be in- 
volved, however, the moment of inertia of the intrinsic 
state is far from sufficient and configuration mixing by 
the diagonalization of Hamiltonian becomes crucial for 
the description of the experimental trend. 



D. Excited states 

We now discuss in Figs. @] and [5] the low-lying states 
for the considered isotopic chains. 

Experimentally (77l - [79| , the excitation energies of the 
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FIG. 2: (Color online) Correlation energy Econ denned in 
Eq. © for Yb, Hf, W, Os and Pt isotopes. 
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FIG. 3: (Color online) Moments of inertia of W isotopes, 
computed by the cranking formula for the coherent state, by 
the the rotor formula L(L + 1) using the 2f eigenenergies of 
the IBM and of the experimental 2+ excitation energies [77I 



ground-state band shown in Fig. 21 namely the 2^, 4+, 
and yrast states, increase as the neutron shell clo- 
sure N — 126 is approached. The increase of these yrast 
levels with neutron number N becomes more rapid with 
smaller Z , when departing from the proton shell closure 
Z = 82. The present results follow the overall experi- 
mental isotopic trend for those nuclei. For Pt, Os and W 
isotopes, the same systematics have been observed with 
the Gogny-DlS functional . 

The LL term has a remarkable influence on the ground- 
state band at the quantitative level. Without this term, 



the experimental yrast spectra would not be reproduced 
with that precision. This is particularly the case with 
lighter W (Hf) isotopes with N = 110 and 112, which 
follow the rotor formula L(L + 1) with their respective 
experimental ratios being E A +/E 2 + —3.27 (3.29) and 
3.23 (3.26) |77j. For these nuclei, the results shown in 
Figs. |4jc) and[4jd) compare rather well with the exper- 
iments. 

We now turn to the description of the side-band ener- 
gies in Fig. [3] To begin with, we discuss the excited + 
( 0^) state. It is well known that the intruder configura- 
tions may play a role for mid-shell Pt isotopes, where the 
oblate-prolate shape coexistence is observed [HI IfiH . The 
phenomenological IBM study (see Ref. [66| , for instance) 
considers particlc-holc excitations across the Z = 82 pro- 
ton shell. In this kind of work one needs to extend the 
boson model space as to take into account the intruder 
configuration with additional proton bosons, arising from 
(mainly) the 2p-2h excitation. The normal and the in- 
truder configurations are mixed, and the model Hamilto- 
nian should be then diagonalized in such enlarged config- 
uration space. The validity of this mixing calculation has 
been discussed extensively [67], 111] , and is thus of great 
interest. 

The mixing in general becomes more significant when 
approaching the middle of the major shell. In Fig. [5ja), 
the calculated 0^ excitation energies for N ^ 116 Pt 
isotopes, as well as those with Gogny-DlS [47|, seem to 
compare reasonably well with the data, even without tak- 
ing into account the mixing between normal and intruder 
states. Furthermore, the original HFB energy surfaces for 
Pt isotopes do not exhibit clear coexisting minima. Due 
to this, the present framework cannot fix the parameters 
for both the normal and the intruder configurations as 
well as those for the operators mixing the two configu- 
rations. Although such a mixing calculation is a rather 
subtle problem, it is very interesting to study the extent 
to which the intruder configuration plays a role when in- 
troduced in the present mapping method. 

It was shown experimentally [loT - [l9j that, in the non- 
yrast states of lighter W, Os and Pt nuclei, the band 
mixing could arise more or less from the coexistence of 
the different intrinsic states mentioned above, and makes 
it rather difficult to identify the clear band structure by 
a model prediction. The band-mixing feature should be 
outside of the model space of bosons with low-spin on 
which the IBM is built, and may be somewhat difficult 
to be reproduced. It is yet not clear whether the similar 
complicated band mixing will be observed in the exotic 
Yb and Hf isotopes. 

The 2^~ level, which is normally the band-head of the 
j^tt _ 2+ (so-called quasi-7) band, is a good test for the 
evolving triaxiality in a given isotopic chain. Figure [5] 
shows that the calculated 2 J level of the N = 116 nuclei 
is lowest among each of Yb, Hf, W and Os isotopes. Ex- 
perimental excitation energies keep steady (decrease) in 
Pt (Os, W) isotopes as N increases from 110 to 116. 

In our calculations, the decrease of the energies of the 
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FIG. 4: (Color online) Theoretical (curves) and experimental [77l - [79| (symbols) low-lying spectra of Yb, Hf, W, Os, and Pt 
isotopes with 110 ^ N ^ 122 for the 2f, 4f , 6^ and 8^ states. Symbols for the experimental levels are defined in the panel 
(c). 



2^~, 3i, 4^ an( i states occurs more rapidly for lower 
Z isotopes, which have a larger number of active bosons. 
Around N = 116 a change in this tendency occurs and 
these excitation energies increase. This is in agreement 
with the only experimental measurement available in Os 
isotopes. 

A remarkable difference between the theoretical and 
the experimental quasi-7-band structure observed in Pt 
and Os isotopes is that the calculated 3^ and the 4^ 
states, and the 5+ and the 6^ states as well, form dou- 
blets, which are absent in the data. Since all the states 
in Figs. EJa) and[5jb) except the 0^" ones, are supposed 
to be the quasi-7 band states, the appearance of these 
doublets points to the emergence of the 7-unstable [80| 
or 0(6) dynamical symmetry (42j . in which the spectra 
belonging to the same family of the quantum number 
r are nearly degenerated. Since the rigid triaxial rotor 
model with 7 = 30° (U predicts the doublets (2+, 3+), 



(4 + , 5 + ), etc, in the 7 band, the experimental data in 
Fig. [5] for (a) Pt, (b) Os, and (c) W isotopes suggest a 
situation rather in between the 7 -unstable rotor and the 



rigid-triaxial rotor pictures. The discrepancy of the 7- 
band energies occurs probably because the IBM energy 
surface does not show the triaxial minimum which is, 
however, seen in the original HFB energy surface. 

There are several possible effects which may eliminate 
this staggering in the 7-band spectra and improve the 
agreement with the experiments at the quantitative level. 
In the present paper, however, we do not look into the 
details of this issue due to the large number of additional 
parameters to be introduced and the lack of experimental 
data for the Yb and Hf nuclei. First, a three-body (cubic) 
term, which partially breaks 0(6) symmetry, may correct 
the deviation. This has been done mainly in the IBM-1 
flkX [75j . For the present case some type of cubic term 
appears to be necessary mainly for W, Os and Pt nuclei, 
where the Gogny HFB energy surface exhibits a shallow, 
but stable triaxial minimum [56j]. While the calculated 
excitation energies of the quasi-7 band for Yb and Hf in 
Fig. [5jd,e) look like that of pure 0(6) limit as well, the 
validity of this term seems to be marginal in these cases. 
Indeed for the Yb and Hf isotopes the original Gogny- 
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DIM energy surface indicates the discrete change of the 
minimum point from the oblate (7 = 60°) to the prolate 
(7^= 0°) sides, similarly to the Gogny-DlS energy surface 

The second possibility would be to relax the con- 
straint on the deformation parameters 7^ and 7„ so that 
they could take different values. As the IBM-2 can be 
viewed as a two-fluid system consisting of proton and 
neutron bosons, the phase-structure analysis would be 
exploited in the context of the coherent-state formalism 
[82| , whereas it is not obvious to define a consistent map- 
ping procedure for realistic cases. 

The third would be the inclusion of higher-spin bosons, 
like the g-boson. It is not independent of the first pos- 
sibility involving the cubic term, since the cubic term 
can be derived effectively from the renormalization of 
the g boson into the sd-boson sector [74j . This would, of 
course, make the problem more complicated. 

We now address the problem of why the side-band 
spectra, particularly for Pt in Fig. [5ja) and Os in 
Fig. Eth) isotopes, are overestimated in the present cal- 
culation when approaching the N — 126 shell closure. 
The direct reason would be that the microscopic Gogny 



energy-surface calculation predicts mostly oblate defor- 
mations with small quadrupole moment but with rather 
large amount of deformation energy characterized by the 
depth of the potential minimum [56| . Such a topology of 
the HFB energy surface is not well described by the IBM 
Hamiltonian close to the end of the major shell Z = 82. 
Nearby the closed shell one has a relatively small num- 
ber of bosons. The deviation of the spectra seems to be 
due to this limited degrees of freedom. The problem on 
the description of the side-band energies was observed in 
other cases of shape transitions in different mass regions 
0, , and is still an open problem. According to the 
above argument it may be expected that the predicted 
levels for exotic Yb and Hf isotopes in the vicinity of the 
shell closure N — 126 might be overestimated. 

To further examine the problem, it is interesting to 
consider the relevant energy ratios, as they nicely trace 
the underlying shape transition. Figure [5] depicts the 
energy ratios (a) i? 4/2 = E(if)/E(2^ ) and (b) i? 47 = 
E(4f)/E(22) as functions of N. The ratio i? 4 / 2 is prob- 
ably the simplest and best-studied measure for the evo- 
lution of collectivity. The ratio i? 47 presents the location 
of the band-head of the quasi-7 band 2+ (2^~) relative to 
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FIG. 6: (Color online) Theoretical (curves) and experimental (symbols) [77| energy ratios (a) R4/2 = E(4f)/ E(2~^) and (b) 
i?4 7 = i?(4f )/I7(2f ) as functions of N. Definitions of the theoretical curves and the symbols for the experimental data appear 
in panel (b). 



the 4 J" excitation energy. Since in many 7-soft nuclei the 
2% level lies quite close to the 4f level, the overall trend 
of ratio i?4 7 can help to measure the 7 softness. 

In Fig. IDJa), the experimental R4/2 ratios for Os and 
W isotopes exhibit a gradual decrease as a function of 
N from the rotor limit of R4/2 = 3.3 in the vicinity of 
N = 110 toward the 0(6) limit of R i/2 = 2.5. This re- 
flects the transition from the axially deformed rotor to 
the 7-unstable shape. Also of particular interest is the 
difference of the R4/2 ratio between Pt isotopes and the 
other isotopes. The experimental R4/2 ratio for all the 
Pt isotopes studied remains practically constant all the 
way, being close to the 0(6) limit of 2.5. The present 
calculation follows the decrease of the experimental R4/2 
value from N = 110 to 116 in Os and W isotopes, while 
an increase is suggested for not only Os and W but Hf 
and Yb isotopes for N ^ 118, contrary to the experimen- 
tal tendency of Os isotopes. The change in the calculated 
ratio R4/2 occurs quite rapidly for Hf and Yb isotopes in 
comparison to W and Os isotopes. The discrepancy of 
the tendency for N ^ 118 for Os nuclei could be the con- 
sequence of the unexpectedly large x-n and \v values with 
positive sign, as seen in Table HI since the corresponding 
IBM energy surfaces exhibit notable oblate deformation. 
The same would hold for explaining the overall devia- 
tion in Pt isotopes. In this context, to describe all the 
observed data including those for TV ^ 118 regime, the 
triaxial dynamics needs to be correctly incorporated in 
the present model. 

The energy ratio i?4 7 is depicted in Fig. HJb). The 
experiment shows that in the lighter Pt, Os and W iso- 
topes with N — 110, 112 and 114, the ratio is below 
unity. While for Pt isotopes the experimental ratio i?4 7 
remains all the way with values close to unity, for Os and 



W isotopes the 7 softness gradually develops with N as 
the ratio i?4 7 increases for 110 ^ N ^ 116 and over- 
passes i?4 7 = latiV = 116. The overall trend of this 
experimental ratio for W and Os isotopes is reproduced 
in the present calculation, and the same systematic trend 
is predicted for Yb and Hf isotopes. For Os, the experi- 
mental i?4 7 ratio decreases from N — 116 to 118, which 
is reproduced by the calculation. In the heavier isotopes 
with N ^ 118 there is a new tendency that the calculated 
ratio shows overall decrease, being much below the unity, 
whereas the experimental ratio for Os isotopes keeps in- 
creasing, being larger than unity. The results presented 
here do not differ much from the case of D1S functional 
already studied in [13, HH . 



E. B(E2) systematics 

Lastly, we examine the B(E2) systematics for a few es- 
sential cases corresponding to the shape transition. The 
B(E2) ratios relevant to the band- head of quasi-7 band, 
2f state, can be the stringent tests. 

We show in Fig. [7] the ratio (a) £(E2;2f 
2f)/5(E2;2 ] f Of) and the branching ratio (b) 

B(E2; 2+ -» 0+)/B(E2; 2+ ->■ 2+) for the considered iso- 
topes in comparison with the data [83l l84j . 

The 2\ — > 2\ E2 transition rate shows a certain sensi- 
tivity to the neutron number N and thus it is useful as a 
signature of the structural evolution involving the 7 soft- 
ness. The B(E2;2+ -> 2+)/£(E2;2+ -> Of) ratios for 
Pt isotopes differ notably from those of other isotopes. 
For Yb, Hf, W, and Os isotopes, the calculated ratio is 
peaked at N — 116. This confirms that in each of these 
isotopic chains the N = 116 nucleus is softest in 7 direc- 
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FIG. 7: (Color online) The B(E2) ratio (a) B(E2; 2+ -> 2+)/B(E2;2+ -> 0+) and the branching ratio (b) B(E2;2+ 
0+)/B(E2;2^ 2f) for relevant low-lying states of the considered Yb, Hf, W, Os, and Pt isotopes with Gogny-DIM EDF. 
Experimental data for W, Os, and Pt isotopes are taken from Ref. [84| . Definitions of symbols and theoretical curves appear 
in panel (b). 



tion. On the other hand, for Pt isotopes the calculated 
B(E2; 2+ 2+)/B(E2; 2+ ->■ Of) value keeps increasing 
toward N = 110 to approach the 0(6) limit, rather than 
taking a maximum at N = 116. This tendency appears 
to be consistent with that expected from the topology 
of the HFB energy surface [55| and from the predicted 
systematics of the quasi-7 band-head in Fig. [5]Ja), which 
reflects that the 7 softness persists for rather wide region 
in the Pt isotopic chain. 

When compared with the D1S case [53j . the present 
DIM result suggests that the ratio B(E2; 2^ -» 
2+)/£?(E2; 2^ — > Of) is rather sensitive to the iso- 
topic chains. In fact, in Fig. [7ta), the £?(E2;2^ — > 
2f)/ J B(E2;2 ] f ->■ Of) values below and above N 
1 1 ti appear to have a certain Z dependence when the 
DIM functional is used. For instance, the -B(E2; 2^ — > 
2^ )/£> (E2; 2+ — > 0+) value for W isotopes is generally far 
from the 0(6) limit all the way. It has been noticed in 
Ref. [53j, however, that the calculated value of this B(E2) 
ratio is practically the same for Os and W isotopes when 
the D1S functional is taken. It would be interesting to 
see if this Z dependence is observed experimentally. 

The branching ratio B(E2; 2+ -> 0+)/B(E2; 2+ -> 2+) 
in Fig. [TJb) also presents a clear signature of the struc- 
tural evolution involving triaxiality. For Yb, Hf and W 
isotopes with 110 ^ N ^ 116, the branching ratio de- 
creases from values close to the SU(3) limit of 0.7 to 
the U(5)/0(6) limit of zero. This behavior corresponds 
to the transition from well deformed to 7-soft nuclei as 
confirmed by the experimental data on Os and W iso- 
topes. At this point, one can observe the increase from 
N = 116 toward the shell closure N = 126. The increase 
represents the deviation from the 7 -soft character, as the 



corresponding mapped energy surface in Fig. [T] exhibits 
notable oblate deformation. The change in the branching 
ratio occurs more slowly than the D1S case [53[. This is 
consistent with our general finding that the DIM energy 
surfaces for these nuclei show less pronounced quadrupole 
correlation than the D1S ones. As observed in Fig. [7f b) , 
the branching ratios for Pt isotopes remain always much 
closer to zero, which is compatible with their sustained 
7-soft character. 



IV. SUMMARY 

In summary, the method of deriving the Hamiltonian of 
the interacting boson model from the constrained HFB 
calculations with the Gogny functional DIM has been 
applied to the spectroscopic analysis of the neutron-rich 
Yb, Hf, W, Os, and Pt isotopes. The microscopic energy 
surface obtained from the constrained HFB calculation 
turns out to be a good starting point for both reproduc- 
ing and predicting the ground-state shape of the consid- 
ered nuclei. Spectroscopic observables that characterize 
the underlying shape transitions, such as excitation en- 
ergies, _B(E2) ratios and correlation energies, have been 
calculated. 

It has been shown that the Pt isotopes largely differ 
from the other isotopes in the rapidity of the shape tran- 
sition. For most of the considered Pt nuclei the mapped 
IBM energy surfaces are 7 soft. The transition occurs 
more rapidly when departing from Z = 76 (Os) through 
Z = 70 (Yb). The triaxial deformation helps to un- 
derstand the prolate-to-oblate shape transition that oc- 
curs in the considered isotopes. The N = 116 nuclei can 
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be commonly identified as the transition points. This is 
most noticeably seen in the overall systematic trend of 
the band- head of the 7 band 2 \ , as well as in energy and 
B(E2) ratios. Predicted spectra have been presented for 
the neutron-rich Yb and Hf isotopes, where a quite rapid 
structural evolution is suggested. When compared to the 
results from the standard Gogny-DlS paramctrization 
|53| . the DIM functional seems to be equally valid 
to describe the physics involved. 

On the other hand, the present work aims at inves- 
tigating the possible ways of refining the current model 
and clarifying its limitations when applied to the consid- 
ered mass region. First, as discussed in Sec. MI Dl the 
discrepancy in the level structure of the quasi-7 band 
turns out to be a major limitation. It is likely that this 
discrepancy is mainly due to the use of the IBM Hamil- 
tonian not reproducing the triaxial energy minimum. A 
specific three-body (cubic) term may improve the agree- 
ment. Second, the boson effective charges need to be 
determined in a microscopic way and effects beyond the 
mean field, like core polarization, should be taken into 
account. It would also be meaningful to compare the 
spectra and the electromagnetic transition rates result- 



ing from the present method directly with those obtained 
from full configuration-mixing and symmetry-conserving 
calculations including triaxial degrees of freedom. This 
would help to quantify the predictive power of the em- 
ployed model when applied to heavy exotic nuclei. Work 
along these directions is in progress. 
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